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Abstract 

We present the analytical results for the K — > 3tt final state interaction phases at 
next-to-leading order (NLO) in CHPT. We also study the recent Cabibbo's proposal 
to measure the tttt scattering lengths combination ao — a<i from the cusp effect in the 
7r°7r° energy spectrum at threshold for K + — > 7r°7r°7r + and Kl — > 7r°7r 7r°, and give 
the relevant formulas to describe it at NLO. We estimate the theoretical uncertainty 
of the ao — 02 determination at NLO in our approach and obtain that it is not 
smaller than 5 % if added quadratically and 7 % if linearly for K + — > 7r°7r°7r + . 
One gets similar theoretical uncertainties if the neutral Kl — > decay data 

below threshold are used instead. For this decay, there are very large theoretical 
uncertainties above threshold due to cancellations and data above threshold cannot 
be used to get the scattering lengths. We do not include isospin corrections apart 
of two-pion phase space factors which are physical. We compare our results for the 
cusp effect with Cabibbo and Isidori's ones. 
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1 Introduction 



Final state interaction (FSI) phases at next-to-leading order (NLO) in Chiral Perturbation 
Theory 1 (CHPT) PS] are needed to obtain the charged CP- violating asymmetries at NLO 
[5j El [7J. The dominant contribution to these K + — > 3tt FSI at NLO are from two-pion 
cuts for topology A in Figured! They were calculated analytically in [5]. 

Though to get the full K — > 3tt amplitudes at order p 6 implies a two-loop calculation, 
one can get the FSI phases at NLO using the optical theorem within CHPT with the 
advantage that one just needs to know 7nr scattering and K — > 3n both at 0(p 4 ). Notice 
that NLO in the dispersive part of the amplitude means one- loop and 0(p 4 ) in CHPT 
while NLO in the absorptive part of the amplitude means two-loops and 0(p 6 ) in CHPT. 

In [5] we took the irn scattering results from [8] and calculated the amplitudes K — > 
3tt at NLO in the isospin limit. We agreed with the NLO K — > 3n results recently 
published in [9] . These were previously calculated in pi)] and used in [TT] , but unfortunately 
the analytical full results were not available there. Now, there is also available the full 
one-loop K — > 3tt amplitudes including isospin breaking effects from quark masses and 
electromagnetic interactions [T2| I13j. 

In the meantime, we discovered some small errata in the published formulas in [5j which 
were corrected in [5J. We correct here another erratum in Eqs. (6.6), (6.7) and (6.9) in [5] 
where we wrote R instead of R. See Section [3] in the present work for the correction. In 
addition, matrix elements 21 and 22 in Eqs. (6.7) and (6.9) in [5] were interchanged. None 
of these errata affects the results for the CP-violating asymmetries nor the conclusions 
given in [3, El [7] . 

The details of the use of the optical theorem to calculate the two-pion cut contributions 
from topology A in Figure [1] were given in Appendix E of the first reference in [5] -solid 
circle and solid square vertices include both tree and one-loop level. We don't repeat these 
details now. 

Notice that diagrams B and C are included in topology A when one takes one of the 
solid vertices at tree level and the other one at one-loop. We plot explicitly diagrams B 
and C since we will refer to them later. 

In the first part of this work, we complete the calculation of all K — > 3ir FSI phases at 
NLO in CHPT by evaluating the three-pion cuts in topologies C, D and E in Figured] for 
K + — > 37r. We also provide with the neutral kaon K^ s ~^ 37r final state interactions at 
0(p 6 ) including both two- and three-pion cut contributions. Thus, we give the full K — > 3ir 
FSI phases at NLO in CHPT in the isospin limit for all K — > 3ir decays. 

The study of FSI in K — > 3n at NLO also became of relevance after the proposal 
by Cabibbo [14J to measure the combination — a2 of tttt scattering lengths using the 
cusp effect in the 7r°7r° spectrum at threshold in K + — > 7r°7r°7r + and Kl — > vr 7r 7r° decay 
ratefl Within this proposal, it has been recently presented in [TB] the effects of FSI at 
NLO using formulas dictated by unitarity and analyticity approximated at second order 

3 Some introductory lectures on CHPT can be found in [T] and recent reviews in [2]. 
4 The cusp effect in SU(2) irir scattering was discussed in [15]. 
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Figure 1: Topologies contributing to FSI in K — > 3n decays at LO and NLO. Squares are 
weak vertices and circles are the strong ones. In topology A, solid circle and solid square 
vertices include both tree- and one-loop level. See text for more details. 



in the irir scattering lengths, aj ~ 0.2. The error was therefore canonically assumed to be 
of order of af, i.e., 5%. There, they used a second order polynomial in the relevant final 
two-pion invariant energy S3 fitted to data to describe the K — ► 3n vertex that enters in 
the formulas of the cusp effect. It is of interest to check this canonical error and provide a 
complementary analysis of the theoretical uncertainty. 

In this work, we use our NLO in CHPT results for the real part of K — > 3ir fitted to 
data to describe the K — > 3ir vertex that enters in formulas of the cusp effect. Notice that 
we do not want to predict the real part of K — > 3tt in CHPT at any order but to have 
the best possible description fitted to data. We treat tttt scattering near threshold as in 
Cabibbo's original proposal. The advantage of using CHPT formulas for the fit to data 
of the real part of K — > 3tt is that it contains the correct singularity structure at NLO 
in CHPT which can be systematically improved by going at higher orders. On the other 
hand, this proposal is not (and cannot be) a CHPT calculation. 

Contributions from next-to-next-to-leading order (NNLO) SU(3) CHPT in the isospin 
limit are expected typically to be around (3 ~ 5)%, so that our NLO results are just a 
first step in order to reduce the theoretical error on the determination of the combination 
a — a 2 to the few per cent level. At NNLO, one can follow a procedure analogous to the 
one we use here to get a more accurate measurement of ao — 0,2 and check the assumed 
NNLO uncertainty. At that point, and in order to reach the few per cent level in the 
theoretical uncertainty, it will be necessary to include full isospin breaking effects at NLO 
too. These are also expected to be of a few per cent as was found for K —> 3n in (T2J [13] . 

In Sections |2] and [3] we introduce the needed notation in the description of K — > 3n 
in CHPT and the K° — ► 3n final state interaction phases at NLO, respectively. This 
completes the K + — > 3n FSI phases already presented in [5]. In Section HI we give the 
relevant formulas of the K — > 3n cusp effect using our NLO formula for the real part of 
K — » 3ir fitted to data. These could be used in a fit of the cusp effect to data in order 
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to extract the ao — a 2 scattering lengths combination as pointed out in [T3] . We compare 
our results with the ones in [16] and give estimates of the theoretical uncertainties in the 
determination of ao — a 2 . We do this both for K + — > 7r°7r°7r + and for Kl — > 7r°7r 7r°. We 
would like to advance here that our analytical results fully agree with those in [16J if we 
make the same approximations. Finally, we summarize our results and main conclusions 
in Section [51 

In Appendix [X] we list the numerical inputs used. In Appendix [B] we present the neutral 
kaon K — ► 3n decay FSI results from two-pion cut contributions. Finally, in Appendix [C] 
we collect the results for all three-pion cut contributions. 



2 Notation 

Here, we want to give the basic notation that we use in the following. In [5] we calculated 
the amplitudes with K m = (k° - /y/2, 
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[^000 ] J 


K 2 (k) - 


-> 7T + (p 1 


)rr-(p 2 )n\p 3 ), 




tfi(fc) - 




)7r~(p 2 )n°(p 3 ), 


W-o 


K+(k) - 


- ApO 


7r°(p 2 )7r + (p 3 ) , 


[^00+] 


K + (k) - 




)-n + {p 2 )^{p 3 ) . 





(1) 

as well as their CP-conjugated decays at NLO (i.e. order p A in this case) in the chiral 
expansion and in the isospin symmetry limit m u = m^. Disregarding CP-violating effects, 
that are negligible in the results reported in this work, K\ ~ K$ and K 2 ~ Kl- 

The lowest order SU(3) x SU(3) chiral Lagrangian describing \AS\ = 1 transitions is 

£\as\=i = C F o e2 °e tr (A 32 utQu) + CT 4 [G 8 tr ( A 32 «X) + G 8 tr (A 32X+ ) 

+ G 27 t ij > kl tr (Ay-u M ) tr (A kl u»)} + h.c. (2) 

with 

C = ~l^V ud V: s ~ -1.07 x 10- 6 GeV- 2 . (3) 

The correspondence with the couplings c 2 and c 3 of [TU] is 

c 2 = CFq Gs; 

c 3 = -^CF 4 G 27 . (4) 
o 

Fo is the chiral limit value of the pion decay constant f n = (92.4 ± 0.4) MeV, 

u M = iu ] (D^U)u ] = u\ , 
Aij = u\ijU ] (Xij) ab = 5 ia 5 jb , 

X+H = wt X wt + (-)ux^u (5) 
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X = dia.g(m u ,md,m s ) is a 3 x 3 matrix collecting the light quark masses, U = u 2 = 
exp (za/2$/Fo) is the exponential representation incorporating the octet of light pseudo- 
scalar mesons in the SU(3) matrix $; 



/ ^ + ^ n + K + \ 

V VE J 



The non-zero components of the SU(3) x SU(3) tensor t 1 -*' are 

1 
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+23,33 



21,13 



= t 
+33,23 



13,21 



3 
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+23,11 _ +11,23 



1 . 
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(6) 



and Q = diag(2/3, —1/3,-1/3) isa3 x 3 matrix which collects the electric charge of the 
three light quark flavors. 

Making use of the Dalitz variables 



x = 



si - s 2 



and 



m 



s 3 — s 



m 



(7) 



with si = (k-pi) 2 , 3s = m 2 K + ml (1) + 
breaking terms] can be written as expansions in powers of x and y 



(2) 



+ m 2 (3) , the amplitudes in (TTJ [without isospin 
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^000 
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+ 0(y 2 ,x), 



(8) 

where the parameters a,, j3i and ji are functions of the pion and kaon masses, Fq, the lowest 
order AS = 1 Lagrangian couplings G$, G' 8 , G27, Ge and the counterterms appearing at 
order p 4 , i.e., L\s and K[s. The definition of these last ones can be found in Section 3.2 of 
[5J. In (jSJ), super- indices R and I mean that either the real or the imaginary part of the 
counterterms appear. 

If we do not consider FSI, the complex parameters c^, /3j and 7$ can be written at NLO 
in terms of the order p 2 and p 4 counterterms and the constants -Bi,o(i) = ^iom + &i 0(1) 
and ^0(1) defined in Appendix B of [5J. There, we gave a f , /3j and 7i at LO in CHPT as 
well as their analytic expressions at NLO. 
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3 FSI Phases for K -> 3 tt at NLO in CHPT 



The strong FSI mix the two final states with isospin 1=1 and leave unmixed the isospin 
1 = 2 state. The mixing in the isospin 1=1 decay amplitudes is taken into account by 
introducing the strong re-scattering 2x2 matrix R. The amplitudes in ([1]) including the 
FSI effects can be written at all orders, in the isospin symmetry limit, as follows [TT] . 
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(9) 



(10) 



projecting the final state with 1 = 1 into the symmetric-non-symmetric basis [TTj. The 
subscript Res (NRes) in ([9]) means that the re-scattering effects have (not) been included. 
In these definitions, the matrix IR, 5 2 and the amplitudes A® depend on si, s 2 and S3. 
Using the usual isospin decomposition of K — > 3tt amplitudes in the isospin limit 



A 00+ (si, s 2 , s 3 ) 
A+^oisx, s 2 , s 3 ) 



^.000 s 2) s 3 ) 
one finds for the elements of 



= 2A c (s u s 2 , s 3 ) + B c (s u s 2 , s 3 ) + B t (si, a 2 , s 3 ) . 
= A c (si, s 2 , s 3 ) - B c (st, s 2 , s 3 ) + B t (si, s 2 , s 3 ) , 

2 r 

= C (si, s 2 , s 3 ) + - 5 4 (s 3 , s 2 , si) - B t (s 3 , si, s 2 ) 

+ A n (si, s 2 , s 3 ) - B n (si, s 2 , s 3 ) , 
= 3A n (si, s 2 , s 3 ) , 
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At LO, R12 = due to the fact that the only contributions to YmAi come from tttt re- 
scattering. At higher orders there are other origins for the contributions to ImAi, therefore 
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Ki2 7^ 0. Notice that the re-scattering matrix R depends on energy and at nir threshold it 
changes -this is in fact the cusp effect discussed in the following sections. 

As explained in [5], we used the optical theorem to calculate the two-pion cut contri- 
butions to FSI phases at NLO for the charged kaon decays which are the dominant ones. 
The analytical result of these contributions for the dispersive part of the amplitudes for 
K + — ► 37r can be found in Appendix E of the first reference in [5]. Here, we also calcu- 
late analytically these contributions to FSI at NLO for all K^ t s ~ > 37r and include the 
three-pion cut contributions in topologies C, D, and E in Figured] in all cases. These last 
three-pion cut contributions have been evaluated numerically -see Appendix O 

The calculation has been done in the isospin limit, apart of the kinematical factors in 
the optical theorem which are taken physical. This constitutes the first full calculation of 
FSI phases in K — > 37r decays at NLO in CHPT. The results for the dispersive part of the 
amplitudes for K^s — > 37r can be found in Appendix [B] for two-pion cuts and in Appendix 
Ofor three-pion cuts for both neutral and charged K — > 3tt. For numerical applications in 
the rest of the paper, we use the inputs in Appendix IA1 

Equations (jHJ) imply the following relation for the two first coefficients of the expansion 
in powers of x and y of the amplitudes in (jSJ 



The matrix R and the phase 62 were given at LO in [5]. We also gave there the two 
combinations of the R matrix elements that can be obtained from the charged kaon decays. 
As we already said in the Introduction, there is an erratum in Eqs. (6.6), (6.7) and (6.9) in 
[5] where we wrote R instead of R. In addition, in Eqs. (6.7) and (6.9) the matrix elements 
21 and 22 were interchanged. 



Recently, Cabibbo showed [H] that the cusp effect in the total energy spectrum of the 
7r°7r° pair in K + — > 7r°7r°7r + is proportional to the scattering lengths combination a® — 0,2 
and proposed to use this effect to measure it. 

This interesting proposal was done at lowest order in the sense that the author just 
considered topologies of type A in Figure [T] with one inr scattering vertex in the final 
state. It has been followed up by a study of higher order re-scattering effects coming from 
topologies B and C in Figure [Tl |16j. The first experimental analysis applying this proposal 
has been already published in p.8J showing clearly that the nominal 5 % theoretical accuracy 
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( 1 + % 82 ) 73,NRes ■ 



(13) 



4 Cabibbo's 



Proposal at NLO 
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quoted in [16] will dominate this determination. It is therefore very important to check 
this 5% theoretical uncertainty and how to reduce it further. 

In [16], in order to make a quantitative evaluation for the Cabibbo's proposal uncer- 
tainties, a power counting in the scattering lengths aj was done. The main conclusion was 
that one needs to include topology C in Figured] to go to 0(a 2 ) accuracy, i.e. around 5%. 
The authors used Feynman diagrammatics to do the counting in powers of dj although they 
did not have an effective field theory supporting it. It is of course interesting to construct 
it and follow that program. In fact, very recently, such an effective field theory has been 
presented [T9] . 

In [TH [16], the real part of the K — > 37r amplitude was approximated by a second 
order polynomial in the relevant final two-pion invariant energy S3 fitted to data. We 
study a variation of Cabibbo's proposal that uses NLO CHPT formulas for the real part 
of K — > 3tt vertex fitted to data instead of the polynomial approximation used in [Lfl [16] 
plus analyticity and unitarity. 

At a given order in the tttt re-scattering, any FSI diagram with at least one two-pion 
cut can be drawn as topology A in Figure [TJ Solid square vertex stands for the effective 
K — > 3tt decay vertex and solid circle vertex stands for the effective tttt scattering. FSI 
diagrams with no two-pion cuts are treated apart. At NLO, these are diagrams D and E in 
Figure [TJ Using analyticity and unitarity, one can cut topology A into two subdiagrams: 
the right-hand side is tttt scattering and the left-hand side is K — > 3tt. Diagrams B and C 
are the only two that appear at NLO in the tttt final state re-scattering, i.e. with two tttt 
scattering vertices and with at least one two-pion cut. 

The cusp effect we are interested in originates in the different contribution of tt + tt~ — > 
tt°tt° scattering to the K + — > 7r°7r°7r + (or Kl — > tt°tt°tt ) amplitude when the invariant 
tt°tt° energy is above or below tt + tt~ production threshold. In Section H~T1 for K + — > tt 7t°7t + 
and in Section H~21 for K° — > 7r 7r°7r°, we obtain these contributions using just analyticity 
and unitarity near threshold, i.e. applying the optical theorem to calculate the imaginary 
part of the discontinuity across the physical tt + tt~ cut and Cutkosky rules to calculate 
the real part of that discontinuity around threshold. In some cases, one can get pieces of 
the real part of that discontinuity by applying the optical theorem below tt + tt~ threshold. 
Sometimes this becomes unphysical because the value of the real world pion masses forbids 
it. In such cases we follow the strategy in [16j of going first to unphysical values of pion 
masses, e.g. m n o > m n +, such that the absorptive part below threshold exists and apply 
the optical theorem in this set up. This result is analytically continued above threshold 
where the amplitude always exists by putting the real value for pion masses. 

As said above, analyticity and unitarity allows to separate tttt scattering, i.e the scatter- 
ing length effects, from the rest in K — > 3tt. In fact, when this tttt scattering is evaluated 
around threshold, one is intuitively lead to use the scattering length to all orders as a 
good approximation for the real part of tttt scattering. Explicitly, we follow [T¥[ IT6] for the 
treatment of tttt scattering matrix elements near threshold. Near threshold, we use [I] 

A 00 : tt°tt° -> 7r°7r°, Re A 00 = 32tt a 00 (s) 
A +0 : tt + tt° -> tt + tt°, Re A +0 = 32vr a +0 (s) 
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A ++ : 
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Re A , 
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(14) 

as a definition of the different effective scattering length combinations a,ij(s), which are the 
unknown quantities to be obtained by fitting the cusp effect. From these fitted effective 
scattering lengths combinations, one can extract the pion scattering lengths by comparing 
them with their CHPT corresponding prediction [16] . Notice that this procedure is sensible 
since the effective scattering lengths are evaluated near threshold and chiral corrections 
are tractable within CHPT. These include radiative and isospin breaking corrections to 
the isospin limit results. These isospin limit results are: a 00 = (a + 0,2)/%, a +0 = a 2 /3, 
dx = (ao — a 2)/3, a + _ = (2a + a 2 )/3 and a ++ = a 2 . We define these isospin limit scattering 
lengths at the following thresholds s t h'- 4m^ + for a 00 , a x , a + _ and a ++ and {ra n + + m^o) 2 
for a +0 . 

In order to have a more accurate description of tttt scattering near threshold than the 
one given by (JHJ), one can follow [16] and perform an expansion around threshold in the 
different kinematical variables on which the amplitudes depend. Up to linear terms, the 
generic matrix elements -neglecting all higher order but the P wave- are of the form [20] 
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4 13 m 2 _ 
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(15) 



(16) 



where ajj(s) are the ones in (114"]) . For numerical applications, we use tq = 1.25 ± 0.04 and 
T2 = 1.81 ± 0.05 [20J which are compatible with the results in [21] and 



ai/2 



(17) 



CHPT 



m 



2 /(12tiFq) to lowest order in CHPT. The rest of af, are zero. 



with a\ 

The above discussion fixes the real part of tttt amplitudes near threshold which are 
the unknowns to be fitted. However, even if S3 is always around threshold, at NLO tttt 
amplitudes are also needed far from threshold and the effective scattering lengths in (fi~4"l) 
do not give a good description of them. As we will see in the next sections, the two cases 
appear clearly separated in the tttt scattering amplitudes at NLO in (1241) and fj4"0]) which 



are either evaluated at S3, at {rn 2 K + 3m^. — s 3 )/2 or at {m 2 K + 7ml + s ^)l^- We only 
leave as unknown the scattering lengths evaluated at S3 around threshold and use the real 
part of the full 7T7t amplitudes at NLO in CHPT expressions in the other cases since the 



5 The effects of tttt threshold singularities are NNLO in a chiral counting and expected to be small. 
However one should include them since as we said we want to treat the tttt scattering non-perturbatively 
- this can be done as in [16] . 
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scattering lengths are evaluated far from threshold. The accuracy in the description of 
the 7T7T amplitudes in these latter cases can be improved systematically by going at higher 
order in CHPT. 

For the other ingredient needed by unitarity, i.e. the real part of K — > 3ir amplitudes, 
the procedure we propose, in order to take into account the singularity structure in the real 
part of K — > 37r amplitudes in a systematically improvable manner, is to treat this real part 
of K — > 37r within CHPT. I.e., first using the tree-level (LO) CHPT formulas, next the one- 
loop (NLO)ones, ■ • •; always fitted to data. This is a difference with the approach in [16] 
where a second order polynomial is used. We want to quantify the numerical differences 
induced by this in the scattering lengths obtained from the cusp. There are also other 
differences coming from some approximations done in [16]. These differences and their 
effects will become clear in the following sections. 

Notice that the use of the CHPT formula for the real part of K — > 37r does not change 
the need of fitting it to data. The aim is not to predict the real part of K — > 3tt but to 
describe its analytic structure as accurately as possible so that, once separated the real part 
of K —>■ 3n from tctt scattering through the optical theorem and Cutkosky rules, one can 
measure nir scattering near threshold as Cabibbo proposed. In particular, tttc scattering 
near threshold is treated non-perturbatively as in [TJ] . The advantage now is that we have 
the correct structure of singularities at a given CHPT order for the real part of K — > 3tt. 
Below we give the formulas which describe the cusp effect at NLO within this approach. 

4.1 Cabibbo's Proposal for Charged Kaon Decays 

Near 7r + 7r~ threshold, we can decompose the K + — > 7r°7r°7r + amplitude as follows [T4"l 116] 



A 



00+ 



A 00+ + Boo+v±(s 3 ), for s 3 > 4m^ + 
A 00+ + iB 00+ v±(s 3 ), for s 3 < 4m^+, 



where Aqo + and -Boo+ are in general singular functions except near the n + ii threshold [22] 
and 

/ \ / |g- (m^ w +m 7r(j) ) 2 | 
Vijis) = d . (19) 

Notice that these kinematical factors are taken with physical pion masses, in this way one 
can describe the cusp effect which is generated by the different behavior for the final two 
neutral pions invariant energy above and below the S3 = 4m^ + threshold. 

With these definitions, the differential decay rate for this amplitude can be written as 

M 

\A 00+ \ 2 = ReA 2 00+ + A A + v±(s 3 )A cusp , (20) 



with 



—r2 



A A = lmA 00+ +v 2 ± (s 3 ) ReB 00+ + lmB 00+ , (21) 
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4 cusp 



-2ReA 00 +Im.Boo+ + 2ImA o+ReS o+ , for s 3 < 4m^+; 
2ReA o+Re-B o+ + 2ImA o+ImS o+ , for s 3 > 4m^ + . 



(22) 



The combination of real and imaginary amplitudes A cusp defined above parametrizes the 
cusp effect due to the 7r + 7r~ — > 7r°7r° re-scattering in the K + — > 7r°7r°7r + decay rate. 

The real part of the amplitude K + — > 7r°7r°7r + , i.e. Re Aoo+, was calculated at NLO 
in CHPT in the isospin limit in [5j [9]. As said before, we don't want to predict it (or any 
NLO counterterm) but to use Re A 00+ at NLO in CHPT formula fitted to data in A cusp . 
This provides a fit at least as good as the second order polynomial parametrization used in 
the original proposal [HI [16], and its precision is just limited by data and the singularity 
structure at NLO in CHPT. 

We get Iim4 o+ and Im_B o+ using the optical theorem, where we consistently use the 
NLO expression for ReA o+ previously fitted to data. 

The singular part of the dispersive amplitude, i.e. ReB 0Q+ v±(s 3 ), is proportional near 
threshold to one half of the discontinuity across the physical cut U which is also fixed just 
by unitarity and analyticity. 

With those results at hand, we can study the proposal of determining the scattering 
lengths combination ao — a 2 from ( 1221 . 

The real part of tctt scattering in the final state interactions is included as in [TS] , using 
just unitarity and analyticity and its treatment was explained in the introduction of Section 
HI In particular, for the real part of 7r + 7r~ — > 7r°7r° scattering near threshold, we use the 
non-perturbative definition a x (s 3 ) in (fT4T) . 

First, we give the expressions for A 00+ and B 00+ using the tree- level CHPT expression 
to fit the real part of K + — > 7r + 7r 7r° to data. We get 



ReA 00+ \ LO (s 3 ) 



ImAoo+Uo(si, s 2 , s 3 ) 



with 



and 



SW^(s) 



PWt°(s) 



+0 

ReB o + \ L o(s 3 
lmB 00+ \ LO (s 3 



+ 



ReG 8 {ml - s 3 ) + 



G 



27 



6(m 2 K 



mi 



[s3(4m^ — 19m 



+5m^ — Ami + ^ m K m ^\ 
voo{s3)a 00 {s 3 ) Rev4 00 + 1 lo(s 3 ) 
v +0 ( Sl )(SW^( Sl ) + (s 2 - s 3 )PW^( Sl )) + Sl ^s 2 
' - - 2ml 



c 



647T/2 
G27 



6(m 2 



K 



mi 



ReG 8 (m 2 K + ml 



-(977-4 + m 2 K (15ml - 19s) + 4m^(s - ml)) 



(s - Am 2 n 

1927T/2 

0, 

C'a x {s 3 



ReG 8 + 



G 



27 



6(m 2 K — ml) 



(19m 2 K - Ami 



1 



(23) 



-ReG 8 (s 3 - ml + m\) + G 27 - [l3s 3 - 7m 2 K - 13m 2 ] 



3 We thank Jiirg Gasser for clarifying this point to us. 
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where C = iC Fq /(/J /#) and the constant C is defined in fl3]). Notice that LO here does 
not mean a calculation at LO in CHPT. It means that we have consistently used the LO 
CHPT expression to fit ReA 0Q+ , while ImA o+ and Im_B 00+ are obtained also at LO in irir 
scattering using the optical theorem. The SW and PW pieces above are the S wave and 
P wave contributions to Iniv4oo+, respectively. 

At NLO in tttt scattering and using the CHPT expression to one-loop [5j [9] to fit the 
real part of K + — ► 7r°7r°7r + to data, we get 



ReA 00+ \ NLO (s 1 ,s 2 , s 3 ) 
lmAoo+\ N Lo(sx, s 2 , S3) 



Rey4 00+ | iO + 



M 7 (s 3 ) + M 8 ( Sl ) + M 8 (s 2 ) + M 9 (s x )(s 2 - s 3 ) 



+M 9 (s 2 )(si - s 3 )] 0(p 4) + SReA 00+ (s 1 , s 2 , s 3 ) , 



ImAoo+lio + woo(s 3 )a o(s 3 ) 



M 7 (s 3 ) + M 8 (s 3 ) 



+M 9 (s 3 )(m^ + 3m 2 n - 2s 3 ) -M 9 s (s 3 



C(p 4 ) 



+ 



ReB 00+ \ NLO (s 3 ) 



v +0 ( Sl )(SW» LO ( Sl ) + (s 2 - s 3 )PW? LO ( Sl )) + Sl ^S 2 
-a 00 (s 3 ) v 00 (s 3 ) lmB 00+ (s 3 ) 

-a x (s 3 ) 2 v 00 (s 3 ) J dxReA 00+ (t + (s 3 ,x),t~ (s 3 ,x), s 3 ) 



-<Ms 3 J - 



dx Iim4 ++ _(t + (s 3 , x), t (s 3 ,x)) 



lm.B 00+ \ NLO (s 3 ) 



Imfioo+Uo + a x (s 3 ) 



2M 11 (s 3 ) + M 10 (s 3 ) + M u (s 3 ) 



-M 12 (s 3 )(m 2 K + 3ml ~ 2s 3 ) + Mf 2 (s 3 ) 



(24) 



when s 3 is near 4m^ + . Where t ± (s,x) = a(s) ± 6(s)x with a(s) and 6(s) defined in 
We have not included the three-pion cut graph contributions which we have checked (see 
Appendix [C] for the results) to be negligible. 

The S and P wave contributions to ImA 00+ -SW+^°(si) and PW+q°(si)- were calcu- 
lated in [5]. With the substitution a (si) — > v + o(si) in the formulas of [5] understood, one 
has 



v +0 ( Sl )SW_ 



^°(si) + S! ^ S 2 
NLO I 



formulas E.25 + E.29 of [4] , 
v +0 (si)(s 2 - s 3 )PW^ LU (s 1 ) + si «-> s 2 = formulas E.26 + E.30 of [4] . (25) 

The functions Mj(s) used in (124")) are the same ones defined in [9j at order p 4 with the 
exception of B(m^^m^^ s 3 ) which has to be exchanged by 



B(m n , m^, s 3 ) = J(m^, s 3 ) 



16tt 2 



log ( ^ ) + 1 



(26) 



with 



1_ f 2 + ^(s 3 )log(i^), for s 3 >Aml + - 



167T 2 1 2 + 2t>(s 3 ) arctan (v(s 3 )) , for S3 < 4m 
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and v(ss) = \/\ss — 4m^ |/s 3 . Functions M»(s) and M/(s) are defined in (|4"T|) as integrations 
on x of M(t + (s, x)) and t + (s, x) M(t + (s, x)), respectively. 

Some pieces of Re B +00 can be obtained going below 7r + 7r~ threshold and applying the 
optical theorem there but one cannot get all of them in such way. The first two lines in 
ReB +0 o i n f|24l) come from diagram B in Figure [1] while the third line comes from diagram 
C in the same figure. 

The real part of the discontinuity across the physical cut in S3 for K + — > 7r°7r°7r + can 
be written as an integral in the x-complex plane between x — — 1 and x = +1 [23J 0. This 
can be expressed as the sum of an integral in x along the real axis with — 1 < x < 1 which 
is finite plus an additional piece which is non-zero for values of S3 above 

'L = ~ m :° m « + -ml + ) (27) 

and diverges as S3 approaches (m^+ — m 7T +) 2 . This piece takes into account the effect of 
the presence of the singularity at the pseudo-threshold s 3 = {m K + — m 7T +) 2 (23] and gives 
an additional contribution to the result for ReB +00 (s 3 ) in (1241) when s 3 is above sl- This 
additional piece is known and can be expressed as the result of an integral over a circuit 
in the complex-x plane around a branch cut. Since we need to describe the cusp just near 
threshold, this extra piece is not needed in ReB +0 o(s3) and one can effectively use the 
result for ReB + oo(s3) in ( 1241) which is therefore finite. This agrees with the naive result of 
applying Cutkosky rules when s 3 below s^. 

We get the function Iim4 ++ _(t 1 , t 2 ) that appears in ReB +00 (s 3 ) from the imaginary 
part of the K + — > 7t + 7t + tt~ decay amplitude, which was obtained in [5] using the optical 
theorem. From that result and disregarding the tiny P wave contribution, we get 

lvaA ++ -(t u t 2 ) = a x {ti)voo(ti) / dy ReA 00+ (t + (ti, y), t~(t u y), h) 



J -1 

+2a+_(t 1 )v±(t 1 ) J dyRel+^it+ih, y),t u t~(ti, y)) 

+ a ++ (t 2 )v±(t 2 ) J ^dyReA ++ _(t + (t 2 ,y),t-(t 2 ,y),t 2 ), (28) 

where the expression for ReA +H at NLO in CHPT can be found in [9] . This expression 

has to be fitted to K + — > 7r + 7r + 7r _ data. Notice that, in the formula above, hit amplitudes 
appear in some cases evaluated far from threshold even if S3 is around threshold. As said 
before, in Section HI these cases are clearly separated and whenever this happens we use 
full NLO CHPT predictions and not the effective scattering length combinations in (|14p 
which are used as unknowns just near threshold. 

The contribution 5ReA 00+ to ReA 00+ \ NLO in the first line of (1241) comes from the dis- 
continuity across the physical cuts and takes into account the singularities of ReAoo + at 



We thank Jiirg Gasser for bringing this work to our attention. 
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Si = (m^ + m n o) 2 (i = 1, 2) and S3 = 4m^. thresholds which start at order p 6 in CHPT. 
Using Cutkosky rules, we get 



5ReA m+ (s h s 2 , s 3 ) = - a o(s 3 )v o(s 3 ) J dx a +0 (t + (s 3 ,x))v +0 (t + (s 3 ,x)) 

x J dyReA Q o + (t + (s 3 ,x),t + (t + (s 3 ,x),y),r(t + (s 3 ,x),y)) 

a +0 (si)w +0 (si) J dx [j 7 00+ {s 1 ,t + (s 1 ,x),t~(s 1 ,x)) 
+v±(r(s 1 ,x))ImB 00+ (r(s 1 ,x))] +S!^s 2 } . (29) 

•7"oo+( s i> s 2, S3) is equal to Iniy4 o+(si, s 2 , s 3 ) in f[211) minus the terms proportional to v +0 (si) 
which produce a regular piece. Again, we know from the result in [23J that (1291) is just 
correct when S3 is in the range 

4m^ <s 3 < — {m 2 K+ - ml + ) (30) 

m^o + 171^+ 

and both S\ and s 2 are in the range, 

(%» + m n +) 2 < si >2 < - (m 2 K+ - m n +(2m n o - m ir +)) . 

(31) 

For larger values of s^, the same comment around (|2TI) applies, i.e., known additional pieces 
that diverge as approaches {m K + —m^i)) 2 have to be added [23]. The discontinuities in 
si and s 2 included in 5Rev4 o+ are completely analogous to the one in s 3 already discussed 
after fl24|) . We do not repeat therefore the discussion already done for the discontinuity in 

S3- 

The difference here is that, while we need the discontinuity in S3 just around its thresh- 
old, it is possible to approach the pseudo-thresholds in si (or in s 2 ) when s 3 is around 
threshold. And when this happens, one has to take into account in SR,eA 0Q+ the additional 
pieces mentioned above which diverge at pseudo-thresholds. A solution to this inconsis- 
tency which does not simply use the discontinuity to describe the cusp effect is discussed 
in [22] . Another possible solution, if one persists in using the discontinuity to describe the 
cusp effect in S3, is to drop out 5ReAoo + and use instead the (unknown yet) full-two loop 
finite relevant pieces to describe the NLO singularities near Sj = {m^+ + m T o) 2 (i = 1,2) 
and S3 = 4m^ thresholds. These additional pieces could also be fitted to data since we 
don't want to predict ReA 00+ but to obtain the best description of it. Here, we don't 
discuss these possible solutions further and leave it for a future study. 

In [TB], they used the approximation that the integrands in the finite integrals (1211) . 
f!25|) and f|29l are to a good accuracy linear in the variables x or y. If we make such 
approximation, we get 

ReB 00+ \ NLO (s 3 ) = ~a 00 (s 3 )v o(s 3 )lmB 00+ (s 3 ) 
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2a x (s 3 ) 2 v 00 (s 3 )ReAoo+ (a(s 3 ), a(s 3 ), s 3 ) 
2a x (s 3 )a x (a(s 3 ))v 00 (a(s 3 ))ReA 00+ (a(a(s 3 )), a(a(s 3 )), a(s 3 )) 
4a rr (s 3 )a + _(a(s 3 ))t; ± (a(s3))Rev4 ++ _(a(a(s3)), a(s 3 ), a(a(s 3 ))) 
2a :E (s3)a ++ (a(s3))f ± (a(s3))Reyl ++ _(a(a(s3)), a(a(s 3 )), a(s 3 )) 



(32) 



and 



5ReA 00 +(si, s 2 , s 3 ) = -4 a 00 (s 3 ) a +0 (a(s 3 )) t» o(s 3 ) w+o(a(s3))ReAoo+(a(s 3 ), a(a(s 3 )), a(a(s 3 ))) 

-2 { a +0 (si)w +0 (si) [j r oo+(si, a(si), a(s x )) + t> ± (a(si)) ImS 00 +(a(si))] 



which agree with [16]. By doing these finite integrals numerically, we have checked that 
the approximation that the integrands are linear in x or y is good to 1 % accuracy. 

Again, in (124)) . the subscript NLO does not mean a calculation at NLO in CHPT. It 
means that we have consistently used the NLO CHPT isospin limit expression to fit the 
effective RqAq Q+ and included im re-scattering vertices also at NLO using just unitarity 
and analyticity. Notice that in particular CHPT is not used to include im re-scattering 
effects. These are unknowns to be determined and are treated to all orders near threshold. 
Sometimes effective irn vertices appear evaluated far from threshold even if S3 is around 
4m^. + . In these cases we use the NLO CHPT prediction and don't leave them as unknowns. 

The relative effect of the cusp in |^4oo+ 1 2 can be seen in Figure [21 where we plot d l rc ^°p( S3 )l 
over using the results in ( 124"1) . Here, 

Tcusp is the contribution of f-i-(s3)A CUS p(si, S3) 
in ([22D to the total K+ -> 7r°7r°7r + decay rate T, 



where N = 5127r 3 m|-, and si max and simm can be found in Eq. (4.4) of [5j. Remember 
that additional pieces to fl29l) and fl33|) discussed after (|3T|) have to be added when either 
Si or s 2 is above the upper limit given in fl3TT) . If added, these pieces would make 5ReA 00+ 
and therefore A cusp (si, s 3 ) to diverge as Si or s 2 approaches {m K + —m n o) 2 . In the integral 
above we do not include those additional pieces which were also not included in [16] and 
leave for a future work the study of the possible solutions mentioned above between fl3"T|) 
and (13"2"1) . So that, we include in A cusp (si, s 3 ) the same terms included in [16J. We perform 
the integrals in x and y in (1241) . fT28|) and (|29l) numerically. 

In Figure [31 we show the cusp effect for K + — > 7r°7r°7r + at NLO. The solid line is our 
result for ^ dr - p ff 3 in (J22D 

using (|24l) . Effective tttt scattering lengths near threshold 
in ( fT4l) are unknowns to be fixed by fitting the cusp effect, however for the numerical 
comparison we take them as inputs given by CHPT [20j at NLO. For the rest of inputs 
needed for ReA 00+ and ReA ++ _ we use the ones in Appendix |AJ which have been obtained 
from a fit to data. 



+Sl 



S2} 



(33) 




(34) 
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d|r cusp (s 3 )|/ds3 
dr( S3 )/ds 3 
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Figure 2: Plot of 100 x d l r ™p( a 3>l over around threshold as a function of s 3 , 4m 2 < 
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|C"| 2 v±(s 3 ) 
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around threshold as a function of S3, 4m 2 < S3 < 4(2m 2 + — 



m^o), for the decay K 
text. 



7r 7T 7T . The meaning of the various lines is explained in the 
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Using the same value of the tttt scattering lengths, we also plot -0™ v±(s 3 ) S3 m Fig ure 
|3]as obtained using the results in [16J -this is the dashed curve. For the slopes needed for 
Rev4 o+ and ReA ++ _ in (4.6) and (4.7) in [16], we use the Taylor expansion of our NLO in 
CHPT results for Rev4 o+ and ReA ++ _. This makes the comparison with the solid curve 
clearer. The effect of using instead the slopes from Taylor expanding |Reyloo+ | 2 + |Imv4 o+ 1 2 
and |Rev4 ++ _| 2 + |Imv4 ++ _| 2 , respectively, amounts to a further increase of the difference 
between solid and dashed curves by around 2 % below threshold and 1 % above threshold. 
Notice that from a fit to data one can just access the slopes from |ReA 00+ | 2 + llmAoo+l 2 
and |ReA ++ _| 2 + |Iim4 ++ __| 2 . 

If one compares the dashed and the solid curves in Figure EJ one gets differences below 
threshold around 3 % while above threshold vary between 2.5 % and 1 %. They come 
from several order 1% approximations done in [16] which we have identified and enumerate 
below. 

The piece Im74 o+ Re5oo+ ; -which we include in the numerics and it was not in [16J- is 
nominally order a 3 and contributes to the cusp below threshold by less than 1 %. Notice 
that this term contains a suppressing velocity factor Voo(s 3 ). Using a quadratic polynomial 
in S3 as in [16] instead of the exact NLO CHPT formula is a very good approximation and 
gives differences smaller than 0.5 %. 

As said above, the difference between doing an integral in x or y in (1241) . (|28|) and 
(|29l) numerically and doing the linear approximation is around 1 % each. Since the cusp 
formula in ( |22l) is quadratic some of these differences are actually doubled. At the end of 
the day, these individually negligible approximations produce the final differences between 
dashed and solid curves quoted above and seen in Figure [31 Notice that we don't get the 
final difference by summing the individual ones, it just happens that they go in the same 
direction. 

This provides a first handle on the theoretical uncertainty of the Cabibbo's approach to 
obtain the scattering lengths since the differences quoted above do not affect the treatment 
of 7T7T scattering. 

Let us now estimate the theoretical uncertainty in our approach of determining the 
effective scattering length a x {s^) in (T14")) . There are two main sources of theoretical uncer- 
tainties. One is how good is the fit of NLO CHPT formulas to experimental data. In case 
of using CHPT NLO formulas for Rev4 o+, this fit produces central values which agree 
with experiment within 3% accuracy [5, 9J. This global fit was done using total decay 
rates and Dalitz variable slopes. We want to remark that this fit of Re v4 o+ to data has 
to be good enough also away from threshold. Notice that, for instance, in (132]) one needs 
ReAoo+ evaluated at a(s-s) which is typically around {m 2 K — m 2 )/2. 

The other main source of uncertainty, which is more difficult to estimate and that is 
present in any approach which describes the cusp effect to NLO, is the NNLO corrections. 
The only definite way to know it is to calculate these corrections. In our approach, this 
means to calculate the real part of K + — > 7r°7r°7r + at two-loop in the isospin conserving 
limit and make the same analysis that we have done but at NNLO. Notice that in this 
analysis the two-pion phase space factors are the physical ones in order to describe the 
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cusp effect. 

Meanwhile, we can just make the following estimate. Going from one order to the next 
one in our approach implies that new topologies with an extra irn scattering vertex are 
needed. For instance, topology C in Figure Q] appears when going from LO to NLO, and 
analogously there appear new topologies when going from NLO to NNLO. Following this 
line, a naive estimate of NNLO re-scattering effects in K — > 3ir is that they are suppressed 
with respect to LO by a?. Notice that the velocity factors that appear after applying the 
unitarity cuts can be order one -see for instance fl32l) . where v±(a(4m*)) ~ 0.6- and do 
not suppress the naive af estimate. Our estimate coincides numerically with the one made 
in [16], i.e. it gives around 5 % corrections from NNLO contributions. As said above, at 
NNLO it is possible to follow a procedure analogous to the one we use here to get a more 
accurate measurement of do — (12 and check the estimated NNLO uncertainty. 

If the theoretical uncertainty from the fit is added to the 5 % of canonical uncertainty 
assigned to NNLO we get a theoretical uncertainty in our approach of obtaining the scat- 
tering lengths from the cusp effect in K + — > 7r°7r°7r + between somewhat larger than 5 % 
(if uncertaintities are added quadratically) and 7 % (if added linearly), i.e., we essentially 
agree therefore with the uncertainty quoted by [IS] . 

A further source of error in the final determination of the scattering length clq — a<i from 
a fit to the cusp effect below threshold is due, as already discussed in [16J, to the existence 
of different strategies to do the fit. The three basic ones are: 

• One can consider all the di £is free parameters in fl24|) . 

• All a,i are fixed to their standard values except the combination — a 2 in a x that is 
extracted from the fit. 

• One can use CHPT as much as possible. Since the combination a —a 2 , or equivalently 
a x( s 3), only appears in -B00+) which is proportional to the cusp, it is enough to keep 
£>oo+ in terms of the scattering lengths. So that, one could fully use the CHPT 
predictions at NLO in [3, [H] for Im A 00+ . 

The comparison of the results for a x obtained from these different fit strategies will provide 
us with another handle to estimate the accuracy of the method. This uncertainty will 
become clear once the fits that we propose are done with data and to be added to the 
previous ones. To the uncertainties discussed above, one still has to add the one here from 
the different data fitting strategies. 

4.2 Cabibbo's Proposal for Neutral Kaons 

In the neutral kaon channel it is also possible to measure the scattering lengths combination 
a — 0,2 from the cusp effect in the energy spectrum of two 7r° of Kl — » 7r°7r 7r°. In this 
case the Bose symmetry of the three neutral pions implies that the amplitude is completely 
symmetric for the interchanges si <-> s 2 «-> S3. The amplitude near 7r + 7r~ threshold can be 
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written as [16 



A 



E [A)oo(-Si) + v±(s i )B 000 (s 1 

i=l,2,3 



if all Sj > Am 2 



000 



E A)0()( s i 
i=l,2,3 



+ zw ± (s fc ) J B oo(sfc) 



(35) 



+ E ^±(si)-B 00 o(si) , 

i=l,2,3 
i==fc 



for s fc < 4m^+ , 
fc = 1,2, 3 

for s^k > 4m^ 



The crucial observation in (1351) is that if the value of Sk is below threshold, the other 
two variables Sj (i ^ k) are of order {m 2 K — m 2 )/2, so safely above threshold. Thus it is 
kinematically impossible to cross the threshold with all three variables at the same time. 
In the region where S3 is around threshold, one can define 



A)ooOi; s 2,s 3 ) 
and #000(33) 



X] ^ooo(si) + ^2 v±(s i )B Q00 (s i 

i=l,2,3 i=l,2 
-£>OOo( s 3) • 



(36) 



Equations (l20l) - (!22l) are now valid also for the decay of K L — > 7r°7r 7r° once the substitutions 

A 00 + -> Aqoo, -4oo+ -> -4ooo and 5 00+ -> # 000 are done. 

Again, we write these amplitudes using the approximation (I14p for the amplitudes near 
threshold. The amplitudes contain 7T7t re-scattering in all channels and in those channels 
where they are far from threshold we use full NLO CHPT expressions which, as already 
pointed out in Section HI is a better approximation than the effective scattering lengths. 

The calculation of v4 000 and -B oo is completely analogous to the charged kaon case - 
see Section H~T1 Using the tree- level result for the real part of Ki — ► 7r 7r°7r°, one finds 



Re^oooUo 
lmA 000 \ LO (s 1 ,s 2 ,s 3 ) 



C (ReG 8 -G 27 )m 2 K , 



ReA 



000 I LO 



[s 3 )a 00 (s 3 ) + ^ v oo(si 



m„ 



i=l,2 



32n 2 f 2 



^ levr/^ 



i=l,2 



R e B 00Q \ LO (s 3 ) 
Im5 oooUo(s 3 ) 



0, 

2a x (s 3 )/ooo( s 3) 



(37) 
(38) 



with 



/ooo(s) 



a 



ReG 8 (s-ml) + 



G 



27 



6(m 2 K 



mi 



(s(9m 2 K - 24ml) ~ 5< + 



+24ml + mlm 2 K ) 



(39) 
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The meaning of LO here is the same that in (|23|) . 

The effect of the charged pion re-scattering appears also at NLO. At this order we use 
the CHPT one-loop formula for ReA 00Q fitted to data. We get 



ReA 000 \ NL o(si, s 2 , s 3 ) 
Im ^oookio(si, s 2 , s 3 ) 



C(P 4 ) 



ReA 000 | LO + [Mo(si) + M (s 2 ) + M (s 3 
+5ReA 000 (si,s 2 ,s 3 ) , 

Imv4 000 | LO + v o(s 3 )a 00 (s 3 ) M (s 3 ) + M (s 3 



+ 



+ 



+ 



ImAjJ' 1 ^!) + W£' 2) ( Sl ) + In^ 6 ' 1 )^) + Im4 6 ' 2 )( Sl ) + 



1(6,2) 



R e -£W) I NLo{s 3 ) 

lmB 000 \ NLO (s 3 ) 



+Sl <-> S 2 

^boo(s 3 ) , 
lmB' 000 \ LO + 2a x (s 3 ) 



M 1 (s 3 ) + M 2 (s 3 ) 



+M 3 (s 3 )(m 2 K + 3ml - 2s 3 ) - M°(s 3 ) 



o(p 4 ) 



where we neglected the contributions of three-pion cut graphs which have been shown to 
be very small in Appendix O The function F qo(s) has the expression 



^000 (s) 



-a 00 (s)v o(s) InxB 000 (s) 
-2a x (s) 2 v 00 (s) |ReA oo (s) + 



dxRe^ooo {t + (s,x)) 



(40) 



— a x (s) / dxlmA + _ (t (s,x),t (s,x)) 



The first two lines here come from diagram B in Figure Q] while the last line comes solely 
from diagram C in the same figure. 

As already discussed in the case of the analogous calculation of the discontinuity for 
the charged K + — > 7r°7r°7r + decay in Section H~T| the integration contour along the real axis 
— 1 < x < 1 in -F 000 (s) gives the correct result just for small values of s, namely [23J, for 



K« < s < - ( 



m 



K° 



m; 



!•)■ 



(41) 



For larger values of s there is an additional piece [23] and the same comments and discussion 
around (T2~Tj) apply here. 

The functions lmA^' l \s) are given in Appendix IB.2I and the meaning of NLO in (j4"U|) 
is the same that in fT24"j) . The imaginary part of the Kl — * K + (pi)7T~(p2)iT (p 3 ) decay 
amplitude, ImA+_ Q (ti, t 2 , t 3 ), has been obtained here using the optical theorem. See Sec- 
tion [3] and Appendix IB.3I for the full expression. From that result, disregarding the tiny 
contribution from the P wave and the pieces that do not contribute to singularities, the 
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function ImA^_ (ti, £ 2 ) is 



ImiJ_ (ti,t 2 ) = a +0 (ti)v + o(ti) / d 3 /ReA + _ (* 1 ,t + (ti,j/),r(ti,y))+ *i -»• h , (42) 



where the expression for ReA^_ at NLO in CHPT can be found in [5], [9] . This expression 
has to be fitted to Kl — > 7r + 7r~7r° data. Notice that, in the formula above, tttt amplitudes 
are sometimes evaluated far from threshold. As explained at the introduction of Section HI 
whenever this happens we use full NLO CHPT predictions for the tttc scattering amplitudes 
and not the effective scattering length combinations near threshold in (|14p . 

The contribution SReA 00Q comes from the discontinuity across the physical cuts and 
takes into account just the singularities of ReA 000 at = 4m^ + , s, = (m n + + m^o) 2 and 
Si = 4m^ (i = 1,2) thresholds. These singularities start at order p e in CHPT. We get 

5ReAo 00 (si,s 2 ,s 3 ) = ^2v±(si) Foooisi) - ^2 v 00 (si) a 00 (si) 



8=1,2 



8=1,2,3 



x / dx < a 00 (t + (si,x))voo(t + (si,x)) ReA 000 (t + (si, x)) 



+ / dyReA 00 a(t + (t + (si,x),y)) + a x {t + (s i} x))v±{t + (s i} x)) 



x dyReA + ^{t + {t + {s i ,x),y),r{t + {s i ,x),y),t + {s i ,x)) }(A3) 



Here, the same comments after (13T!) apply. I.e., the integration contour — 1 < x < 1 gives 
the correct result just when both S3, s\ and s 2 are all in the range in (T4T!) -remember that 
S3 is always around threshold. For values of s% or s 2 outside (j4"Tl) the same discussion after 
( IHTj) applies. 

In [TB], they used the approximation that the integrands in (j4"0"j) . (j4"2"j) and (T4"5|) are to 
a good accuracy linear in the variables x and y . If we make such approximation, we get 



R e B Q0Q \ NLO (s 3 ) 



-a o(s 3 )v 00 (s 3 )lmB 000 (s 3 ) 

-2 a x {s 3 ) 2 v 00 {s 3 ) { ReA 000 (s 3 ) + 2 ReA 000 (a(s 3 )) 



5ReA 000 (si,s 2 ,s 3 ) 



-8a :c (s3)a + o(a(s3))w + o(a(s3))ReA + _ (a(s3), a(a(s 3 )), a(a(s 3 ))) 
-JZ w± ( Si ) {aoo(si)woo(si)Im J B 00 o(si) + 

i=l,2 

+2 a x (sj) 2 f o(si) Rev4 000 (s;) + 2 ReA 00 o 

+8a a ,(s i )a + o(a(sj))w + o(a(si))Re^_ (a(s i ), a(a(s;)), a(a(sj)))| 
-2 aoo(s<)voo(s«){aoo(a(si))voo(o(s<)) [^A 000 (si) + 2ReA 

i=l,2,3 

+2a a .(a(s i ))u±(a(s i ))ReA^_ (a (a(s;)) , a (a(sj)) , a(s;)) j 



000 
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lOOx 



d|r cU sp(s3)l/dS3 




over 



dr( S3 ) 



Figure 4: Plot of 100 x 

s 3 < 4(2m 2 + — m 2 ), for the decay K L — > 7r°7r°7r 



0.078 0.08 0.082 S 3 (GeV 2 ) 
around threshold function of S3, 4m 2 < 



which agree with [16] up to the terms proportional to Voo(a(si)) and v±(a(si)) in SReA 000 , 
which were not included there. 

In Figure H] we show the ratio 100 

x d IW*3)l over dT^) uging 

the results in (|4"0~]) for 

K L — > 7r°7r°7r . Here, r cusp is the contribution of t>±(s 3 )A cusp (si, s 3 ) in (1221 to the total 
Kl — > 7T 7r 7r° decay rate T. We perform the integrals in x and y in (|40p . (|42p and (|43p 
numerically. 

For the sake of numerical comparison with the results in [16], we do not include the 
additional pieces discussed after (14T1) in the integral over si in (|34|) defining d|r cusp (s3)|/ds3, 
which were also not included in that reference either. An analogous discussion to the one 
after (1541) applies in this case as well. 

In Figure we plot the cusp effect ^™ ( p ff 3 in J22J for K L -> ttVtt at NLO. The 
solid line are our results in ( 1401) . We use the same inputs as for the charged case in the 
previous section. 

We also plot in Figure El the cusp effect -ML dr ^ 5 ^) S3 using the results in [16] -this 
is the dashed curve. For the slopes needed for ReAQ 0Q and ReA^_ in (4.58) and (4.59) 
in [IE], we use the Taylor expansion of our results for ReAQ 0Q and ReA+_ . This makes 
the comparison with the solid curve clearer. The numerical differences between solid and 
dashed curves increase by around 2 % if instead we use the slopes from Taylor expanding 
I Retool 2 + llmA^ol 2 and |Rev4^_ | 2 + |Imv4^_. | 2 , respectively. 

If one compares the dashed and solid curves below threshold in Figure O one gets 
differences between 3 % and 1 %. As for the charged case, we have identified the origin 
of these differences: they come from several order 1% approximations done in [TB] . Notice 
that we don't add them, it just happens that they go in the same direction. The part of 
the singularities in 5ReA 0Q0 that was not included in [16] -see comment after (145!) - and 
that we do include in the solid line contributes about +1.5 % to this difference. The other 
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3N dT cusp /ds 3 (r* c y6 
\C'\ 2 v ± {s 3 ) V 



0.01 



\ 0.074 0.076 0.p78 0.08 0.082 
\ 



-0.01 



-0.02 



-0.03 




s 3 (GeV 2 ) 

Figure 5: Plot of |§7p ^ff* 3 around threshold as a function of S3, 4m 2 Q < S3 < 
4(2m^ + — m^o) for the decay — > 7r 7r 7r°. The meaning of the various lines is explained 
in the text. 



approximations are analogous to the ones already enumerated for the charged case -see 
previous section. For instance, the piece IulAqqq Rei? 000 -which we include and it was not 
in [IS]- is nominally order a 3 v±(s3) and contributes to the cusp below threshold by around 
1 %. Again, at the end of the day, these individually negligible approximations produce 
the final differences between dashed and solid curves quoted above. 

Notice again that the differences between our approach and the one in [16] do not affect 
the treatment of the nir scattering part. 

Above threshold, there are large numerical cancellations between Imy4 000 Im5 000 and 
Rey4 000 Re_B 000 in the corresponding expression for A cusp in Kl — > 7r 7r 7r° -see (1221) for the 
analogous charged case. As a consequence, one should know both the real and imaginary 
parts of A 0Q0 and B 000 with a precision better than 1 % to predict the cusp above threshold 
with an uncertainty around 5 %. The observed difference between the two curves above 
threshold in Figure is due to the need of this fine tuning and is not physical. This 
fact makes the region above threshold not suitable to extract the scattering lengths in 
K L -> 7r 7r°7r°. 

Let us now estimate the theoretical uncertainty in our approach of determining the 
effective scattering length 0^(53) in (IT41) from Kl — ► 7r 7r°7r°. As for the charged case, one 
has the theoretical uncertainty from the fitting to ReA oo which includes the theoretical 
error which measures the accuracy of the formula used to do the fit. As said above, this 
has to be checked once the real fit is done but we believe that it is realistic to assume that 
this theoretical uncertainty is around 2 %. Again, the accuracy of the data should be at 
the level of a few per cent as for the charged case. 

If this is added to the 5 % of canonical uncertainty assigned to NNLO we get a theo- 
retical uncertainty for the scattering lengths combination a — 0,2 from the cusp effect in 
Kl — > 7r 7r 7r between somewhat larger than 5 % (if uncertainties are added quadratically) 
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and 7 % (if added linearly). I.e., we predict a similar theoretical accuracy in the extraction 
of 

a o ~ a 2 from neutral and charged kaon cusp effects. But notice that in the neutral case, 
this uncertainty just applies to the analysis of the data below threshold. As said before, 
there are very large numerical cancellations above threshold which preclude the use of these 
data. 

To the uncertainties discussed above, one still has to add the one from the different 
data fitting strategies as described in Section 14.11 

5 Summary and Conclusions 

In Section [3l we have presented the full FSI phases for all K — > 3n decays at NLO in CHPT, 
i.e. at 0(p 6 ) analytically. The two-pion cut contributions for K + — > 3n were already 
presented in [5]. We complete the calculation here with the three-pion cut contributions 
and the full result for Kl,s ~^ 37r. The two-pion cut contributions are given analytically 
while the three-pion cut ones are done numerically and checked always to be negligible. 
We used the techniques already explained in [5] which are based on perturbative unitarity 
and analyticity of CHPT. 

In Section 14.11 we study Cabibbo's proposal to measure the scattering lengths combi- 
nation ao — ct2 from the cusp effect in the total 7r°7r° pair energy spectrum in K + — > 7r 7r°7r + 
[HI [16]. To be more specific, ours is a variation of the original Cabibbo's proposal that 
uses NLO CHPT for the real part of K — > 37r vertex instead of the quadratic polynomial 
in S3 approximation used in [EJ [16] plus analyticity and unitarity. We studied also the 
analogous proposal for K L — > 7r 7r°7r° in Section 14.21 

Notice that we do not use CHPT to predict the real part of K — > 3rr, but use its exact 
singularity form at NLO in CHPT to fit it to data above threshold. If the two-loop CHPT 
singularity structure was known it could be used in order to take this singularity structure 
exactly in ReA o+- The treatment of 7T7t scattering near threshold is independent of this 
choice and we treat it in the same way as in |16| . See the introduction of Section [H 

The cusp effect originates in the different contributions to K + — > 7r°7r°7r + and — > 
7r 7r°7r° amplitudes above and below threshold of 7r + 7r~ production in the 7r°7r° pair invariant 
energy. We obtain these contributions using just analyticity and unitarity, in particular 
applying Cutkosky rules and the optical theorem above and below threshold to calculate 
the discontinuity across the physical cut. This allows us to separate tttt scattering -which 
we want to measure- from the rest of K + — > 7r°7r°7r + or Kl — > 7r 7r 7r°. 

We would like to remark here that making the same approximations that were done in 
[T6] we fully agree with their analytical results. In particular, we checked that the use of 
the quadratic polynomial in S3 in [16] instead of CHPT formulas at NLO produce negligible 
differences -around 0.5 %- in A cusp in (122]) . 

The real part of the discontinuity has a singularity when any of the Sj invariant energy 
reaches its pseudo-threshold at {m K — m^i)) 2 as described in |23j. We have discussed how 
this singularity appears in our formulas for the discontinuity and discussed its effects in 
the description of the cusp effect using the discontinuity -see Section 14.11 and Section 14.21 
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In particular, we pointed out that while the presence of that singularity at pseudo- 
thresholds does not affect Rei?oo+(s3) in (|24T) and Re-B 000 (s3) in (f40|) when S3 is around 
threshold, one needs to take fully into account its effects for 5ReA 00+ in (|29|) and (1331) and 
5ReA Q0Q in (H3|) and (145!) when si or s 2 is above {rn 2 K — m^)/2. For a possible solution of 
this problem which does not simply use the discontinuity to describe the cusp effect see 
[22] . Another possibility could be to use, instead of 5ReA 00+ and 5ReA Q00 , the full-two 
loop (not available yet) finite relevant pieces to describe the singularities at thresholds at 
NLO in ReA 00+ and Rej4 000 , respectively. This could be fitted to data since we don't want 
to predict Rev4 o+ and Rev4 000 but to obtain the best description of them. We leave a more 
detailed study of this problem and possible solutions for a future work. 

In Sections 14.11 and 14.21 we have also discussed some approximations done in [16] and 
the numerical differences they induce in A cusp in (j221) . See Figure [3] and Figure [5] and text 
around. We have identified them and found that though each one of them is individually 
negligible (between 0.5 % to 1 %) they produce final differences in the A cusp around 3 %. 

In the same sections, the theoretical uncertainties in the determination of clq — a>2 if 
one uses our formulas to fit the experimental data are discussed. We concluded that for 
K + — > 7r°7r°7r + , this uncertainty is somewhat larger than 5 % if uncertainties are added 
quadratically and 7 % if added linearly. I.e., we essentially agree with the estimate in [16]. 
Notice that we get our final theoretical uncertainty as the sum of several order 1% to 2% 
uncertainties to the canonical NNLO 5 % uncertainty. 

For the case Ki — > 7r 7r 7r°, we get that -if one uses just data below threshold- the 
uncertainty in the determination of ao — a 2 is of the same order as for K + — > 7r 7r 7r + . 
Above threshold, we found large numerical cancellations which preclude from using it. 

An expansion in the scattering lengths cij and Feynman diagrams were used in [16] to 
do the power counting and obtain the cusp effect description of K + — > 7r°7r°7r + at NLO. 
In general, when FSI nn scattering effects are included at n-th. ordeJ§, there appear new 
topologies in K — > 3n -topology C at NLO order, for instance- which give contributions 
to A cusp of order a™. The canonical uncertainty of the n-th order results is thus a™. Notice 
that the velocity factors that appear after applying the unitarity cuts can be order one -see 
for instance (1521) . where v±(a(4m^)) ~ 0.6- and do not suppress the naive a™ estimate. 

Our estimate for the uncertainty from NNLO, ~ a], coincides numerically with the one 
made in [IB] .i.e. it gives around 5 %. We conclude that one cannot expect to decrease this 
canonical 5% theoretical uncertainty of the NLO result unless one includes tttt scattering 
effects at NNLO. If one wants to reach the per cent level in the uncertainty of the determi- 
nation of a — a 2 from the cusp effect, one would need to include those NNLO re-scattering 
effects. As said above, at NNLO it is possible to follow a procedure analogous to the one 
we use here to get a more accurate measurement of ao — a 2 and check the estimated NNLO 
uncertainty. 

We have just included isospin breaking due to the different thresholds using two-pion 
physical phase spaces in the optical theorem and Cutkosky rules. This is needed to describe 
the cusp effect. The rest of NLO isospin breaking is expected to be important just at NNLO. 

8 n = 1 order stands for LO contributions. 
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At that order, isospin breaking effects in tttt scattering at threshold -both from quark 
masses and from electromagnetism- will have to be implemented and their uncertainties 
added. 

Finally, we believe that it is interesting to continue investigating the proposal in [TU [16] 
to measure the non-perturbative nn scattering lengths from the cusp effect in K + — > 
7r 7r°7r + and Kl — > tt tt tt . Another interesting direction is to develop an effective field 
theory in the scattering lengths which could both check the results in [16] and allow to 
go to NNLO. This type of studies is already underway and firsts results were presented 

P2H22]. 
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Appendices 

A Numerical Inputs 

Here we discuss the numerical inputs we use in the numerical applications. In [9], a fit to 
all available K — > tctc amplitudes at NLO in CHPT and K — > 3tt amplitudes and slopes 
in K — > 3tt amplitudes was done. This was done in the isospin limit. More recently, 
this fit was updated in [T3] with new data on slopes and including also the full isospin 
breaking effects. Though our calculation of FSI at NLO uses isospin limit results we will 
use the results obtained in this last fit. The reason is that the change in the fit results is 
due to both the new data used and the isospin breaking corrections at the same level and 
therefore cannot be disentangled. In addition, the main isospin breaking effects due to the 
kinematical factors is also taken into account in our results. 

At order p 4 , NLO in CHPT, the results in [13] are equivalent [using F = 87.7 MeV] to 

ReG 8 = 6.6±l.l and G 27 = 0.44 ± 0.09 . (46) 

In this normalization, ReGg = 1 = G27 at large N c . (N c is the number of colors of QCD). 

For the NLO prediction of the FSI, we only need the real part of the counterterms and 
in particular the combinations in Table 1 of [5] which from the new fit in [13] are given 
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ReKi(Mp) from [13] 


K 2 (M P ) 


G 8 x (48.5 ± 2.4) • 10~ 3 


K 3 (M p ) 


G 8 x (2.6 ± 1.2) ■ 10- 3 


K 5 (M P ) 


-G 27 x (41.2 ± 16.9) • 1(T 3 


K e (M p ) 


-G 27 x (102 ± 105) • 10" 3 


K 7 (M P ) 


G 27 x (78.6 ± 33) • 10~ 3 



Table 1: Results for the order p 4 counterterms KeKi from the fit to data done in [T3] . 
The values of Re Ki which do not appear are zero. For definitions of the counterterms, see 



in Table \M These were obtained from a fit of the NLO in CHPT results to experimental 
data. Since we only want a good fit to data, one can fix them to any of these values which 
produce equally good fits. 

B FSI Phases at NLO for Neutral Kaon Decays: Two- 
Pion Cuts 

Here we include the analytical results for the two-pion cuts contributions to the dispersive 
part of the neutral decay amplitudes at NLO in CHPT, i.e. 0(p 6 ), coming from diagrams 
in Figure [MZ1 Analogous results for the charged decay amplitudes as well as a more 
detailed description of the method, based on the use of the optical theorem, can be found 
in Appendix E of the first reference in [5]. The fully NLO FSI phases are completed by 
the calculation of the three-pion cut contributions, whose analytical results are given in 
Appendix O In Subsections IB.2[ IB.3l and lB.4l we give the dispersive part of the amplitudes 
Ai for Kl s 3vr denoted by Im A\ where used the super-index 6 to indicate the CHPT 
order. 



B.l Notation 

In all the definitions and results written in this Appendix and in Appendices IB.2[ IB. 31 and 
IB.4[ the functions Mi(t) and Pi(t) are only the real part of the corresponding functions 
defined in [9] and [8] respectively. The expressions in Eqs. fl55|) - fl70|) are thus real. 
We define 

Mi{s) = j dx Mi(a(s) + b(s)x)\ p4 , 

M°{s) = J dx (a(s) + b(s)x)Mi(a(s) + b(s)x)\ p4 , 
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Mf s (s) 



dx (a(s) + b(s)x) 2 Mi(a(s) + b(s)x)\ 



(47) 



with a(s) = -(3s - s) H — . 



and b(s) 



1 

2~7 

x 



(s - (m w (i) + m 7r( 2)) 2 )(s - (m^ti) - m W(2) ) 2 ) 



x 



(s - (tok + m 7r(3 )) 2 )(s - (m K - «V 3) ) 2 ) 



1/2 



(4* 



for s = (k — P3) 2 , see (pQ) for definition of these momenta. 

The amplitudes at 0(p A ) for the 7r7r — > 7r7r scattering in a theory with three flavors can 
be found in [8|. We decompose the amplitudes in the various cases as follows. For the case 
7r + 7r + — > 7r + 7r + the amplitude at C(p 4 ) is 



Ui = P 1 (s) + P 2 (s,t) + P 2 (s,u). 
For the case 7r°7r° — > 7r + 7r~ the amplitude at 0(p ) is 

n 2 = p 3 (s) + p i (s,t) + p i ( s ,u). 

For the case 7r + 7r~ — > 7r + 7r~ the amplitude at 0(p ) is 

n 3 = p 5 (s) + P6(s,t) + P 6 (s,M) + p 7 ( S! t)-p 7 ( s , M ). 

Finally, the amplitude 7r 7r — ► 7r°7r° at 0(p 4 ) is 

n 4 = p 8 (s) + p 8 (t) + p 8 (u). 

The value for the various Pi can be obtained from [S|. In the following we use 



(49) 



(50) 



(51) 



(52) 



Pl n )(s) = c(s) n I dx(l-x) n P(c( S )(l+x),c( S )(l-x)), 



s n c(s) m 



dx (l-x) m Pi(s,c{s)(l-x)), 



with c(s) 



B.2 FSI for -> 7r u 7r u 7r u at NLO 



c( S )"y dx{l-x) n Pi(c(s)(l-x) lS ), 
~(4m 2 -s). 

o^o^o 



(53) 
(54) 



Diagrams A -two charged pions in the loop- and B -two neutral pions in the loop- in 
Figure [6] correspond to the two possible contributions for Kl — > 7r°7r 7r°. We first compute 
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71 1 
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B 



Figure 6: Diagrams for the calculation of FSI for Kl — > 7r°7r 7r° using the optical theorem. 
The square is the weak vertex and the circle is the strong one. 

the case when the weak vertex in Figure [6] is of 0(p 4 ) and the strong vertex of 0(p 2 ). The 
results are 

lmA^ 1 \s 1 ,s 2 ,s 3 ) = lmA^ 1 \s 1 )+lmA^ 1 \s 2 )+lmA^ 1 \s 3 ), 

lmA^ 1] (s) = ^^{s - ml) [Mi(a) + M 2 (s) + M 3 (s)(m 2 K + 3m 2 - 2s) - M s 3 {s) 

ImA§ 2 \ Sl ,s 2 ,s 3 ) = hnA$ 2 \ Sl )+lmA% 2 \s 2 )+lmA$ 2 \s 3 ), 

,(6,2)/ v _ O-(s) 2 



327T/2 



0(p*) 



Then we consider the same diagrams of Figure [6] with a weak vertex of 0(p 2 ) and a strong 
vertex of 0(p A ). We have 

Im^K^ss) = lmA^(s 1 )+ImAi 6 ' 1 \s 2 )+lmA^ 1 \s 3 ), 



a(s) 
IQtc 

hnA$>*>(si, 5 2> s 3 ) = Im4 6 ' 2 )( Sl ) + Im^ 6 ' 2 ^) + Im^ 6 ' 2 ^ 



Im4 6 ' 2 )( S ) = ^ (Mi(s) + M 3 (s)(m 2 K + ml - 3s))\ (P 3 (s) + P 4 (0 ' 0) ( S )) 

107T P 



MA^(a) = (M (s 3 ) + M ( Sl ) + M (s 2 ))\ p2 (P 8 (s) + lf^\s)) , (56) 

for those diagrams with two charged and two neutral pions in the loop respectively. 

B.3 FSI for K L -> tt+tt tt° at NLO 

The calculation is analogous to the one for Kl — > 7r 7r 7r°. The three contributions in 
Figure [7] correspond to two charged pions in the loop -results in Eqs. (j57jl and flSTj) . two 
neutral pions in the loop -results in Eqs. f)58j) and fl62l) - and loops with one neutral and 
one charged pions. In the last case we have a S-wave contribution -results in floTJl) and f[6"31 . 
and a P wave contribution -results in fl6"U|) and (1M1) . respectively. Eqs. (JBTjl . (JSED, dSHD, 
and (I60I) are the results of using the weak vertex at 0(p A ) and the strong vertex at 0(p 2 ). 



0(p*) 



ImA^( S ) = -^miMoW + AfoW^,- (55) 
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B 



Figure 7: Diagrams for the calculation of FSI for Ki — > 7r + 7r 7r° using the optical theorem. 
The square is the weak vertex and the circle is the strong one. 



Eqs. (ISTl) . fl62l . (163]) . and fl64|) are the results of using the weak vertex at 0{p 2 ) and the 
strong vertex at 0(p A ). 



ImA 



(6,1) 
W 



ImA 



ImA 



(6,2) 
W 

(6,3) 
W,S 



Imyli 6 ' 1 ) 



7T,5 



167T/2 

-M*(s 3 ) 



(s 3 ) Mi(s 3 ) + M 2 (s 3 ) + M s (s 3 )(m 2 K + 3ml - 2s 3 ) 



J 0(p*) 



327T/2 



(s 3 -m 2 ) M (s 3 ) + M (s 3 ) 



J 0(p*) 



(m 2 - Sl /2) 2M 2 (s 1 ) + M 2 ( Sl ) + M 1 (s 1 ) 



32tt/2 

-M 3 (si)(m^ + 3m 2 - 2si) +M|(si) 



o(p*) 



Si <-> s 2 



(57) 
(58) 

(59) 



o-(si) 



si(s 2 - s 3 ) 



647r/ 2 s 2 - 2si(m^ + m 2 ) + (m^ - m 2 ) 2 
( Sl - - 3m 2 )(M 2 ( Sl ) - Mi(si) 

+(2 Sl -m\- 3m 2 )M 3 ( Sl )) + (2M|( Sl ) - 2M{(s x ) 

+ (5 Sl - 3m| - 9ml)M°( Sl )) + 2M? {s x ) + h( Sl ) 2 M 3 ( Sl ) 

+Si <-> s 2 . 

0"(g3) 

167T 
0"(g3) 

327T 

g ( g l) 

32tt 



o(p 4 ) 



(60) 



(M^sa) + M 3 ( S3 )(m 2 V + 3m 2 - 3s 3 ))L (P 5 (s 3 ) + ^(^s)) , (61) 



(M (s 3 ) + M (si) + M (s 2 ))L a (P 3 (s 3 ) +P 4 (0 ' 0) ( S3 )) 



M^Si) + ^M 3 (s 2 )( Sl - s 3 ) + ^M 3 (s 3 )( Sl - s 2 ) 



x (P 3 (0 ' 0) ( Sl ) + P^(s x ) + P 2 \ u 4 ; ( Sl ) ) + Sl <- s 2 



>(0) 



(0), 



(62) 



(63) 
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Im4 6 ? = ^(3M 3 ( Sl )(5 2 - 5 3 ))| J 



64vr v v yv Jnp Si- Ami 

x - ^(P^i) _ p( )( Sl ) + i5(0) (si)) + 2 p(i,o )(si) 

-2P$( Sl ) + 2P 2 ( i 1 4 ) (s 1 )) + Sl <- s 2 . (64) 



B.4 FSI for K s -> tt+tt tt° at NLO 

The diagrams contributing to this decay are the same depicted in Figure [7J The result 
corresponding to diagram A in Figure [7J is in Eq. [65] for the weak vertex at 0(p 4 ) and the 
strong vertex at 0(p 2 ) and in Eq. [68] for the weak vertex at 0(p 2 ) and the strong vertex 
at 0(p A ). In this case there is only P wave contribution. 

Diagram B in Figure [7J does not contribute. Diagram C in Figure [7J gives both an 
S-wave - results in Eqs. (1661) and (|69|) - and a P wave contribution - results in (!67j) and 
(ITU]) . Equations (|66l) and (l6Tj) are the results of the case in which the weak vertex is at 
0(p 4 ) and strong vertex at 0(p 2 ) in diagram C. Equations ([69]) and (ITU]) are the results for 
the case in which the weak vertex is at 0(p 2 ) and the strong vertex at 0{p i ) in diagram 
C. 

lmA WP = tMz 2 , 2 l 3iS 2\~ S l\ 2 ^|a(s 3 )(M 4 (s 3 ) + (2a(s 3 ) 

lQTif 2 s 3 -2(m 2 K + ml)s 3 + (m 2 K -ml) 2 I 

-s 3 )M 5 (s 3 )) - (MZ(s 3 ) + (3o(s 3 ) - s 3 )M£(s 3 )) + Mf s (s 3 ) 
_^( S3 ) 2 M 6 (s 3 )j (65) 

g ( g i) 

64tt/= 

-M 6 (si)) - M|( Sl ) + M|(si)} - si <- s 2 , (66) 



ImAg = |M( 2m ' - Sl )|2M 4 ( Sl )-M 4 ( Sl ) + (m^ + 3m2 -2 Sl )(M 5 ( Sl ) 



lmA<$ 



32nf 2 sf - 2s 1 (m 2 K + ml) + {m 2 K - ml) 



x{ l&( Sl ) 2 M 5 (si) + a(si)(M 4 (si) - (2a(si) - 8 1 )(M s {s 1 ) + M 6 ( Sl ))) 



M 4 s ( Sl ) + (3a(si) - Sl )(M 5 s ( Sl ) + M 6 s ( Sl )) - Mf( Sl ) - M 6 ss ( Sl ) 
si <-> s 2 (67) 

(s 3 -4m2)Pf )( S 3) + 2p(^)( S3 )" 



T ,(6,1) a(s 3 ) (M 4 (si)-M 4 (s 2 ))| p2 

llll/ 1 d — 7i 

n ' p 16tt s 3 - Ami 



(68) 

= ^M 5 (s 1 )\ p 4Ss 1 -(m 2 K + 3m 2 M^ '°\s 1 ) + P^ 

+P[%i))-si^ s 2 , (69) 
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Figure 8: Example of diagram for the calculation of FSI for K + — > 7r°7r°7r + with momentum 
assignment. 

x (( Sl - Aml)(P^\ Sl ) - P$( ai ) + Pi°2( Sl )) + 2P 3 ^)( Sl ) 

-2P$( Sl ) + 2P 2 ( ) 1 2(si)) - «i <- a 2 ■ (70) 

C Three-Pion Cut Contributions to FSI Phases 

To calculate the contributions of the topologies in Figure [TJ C, D, and E to the FSI, we 
need the tree level vertices of K — > 3ir and 3tt — > 37r. The complete tree level amplitude of 
the 37r — > 3n scattering can be easily calculated using the code Ampcalculator developed 
in [21]. In this way all possible 3n — > 3n are included. In order to perform the integral in 
the phase space as prescribed in the optical theorem we assign the momentum to pions as 
shown in an example in Figure [H The momentum assignment is done preserving the suffix 
3 for the odd pion in the K + — > 3tt vertex. 

In order to perform the integral in the phase space over momentum it is necessary 
the fix a reference frame. We chose the reference frame in which the decaying kaon is at 
rest and the momentum p 3 defines the z-axis, 

p K = (M K , 0,0,0), P3 = (p§,0,0,p 3 ), 
Pi = (Pi, 0, pi sin rj,p 1 cosrj), p 2 = -Pk -Pi -Pa • (71) 

The momenta describe also a decay of a kaon into three pions. The most general 
expression for these momenta, in the chosen reference frame, is 

qi = R z (a)R y (P)R z (i) n , 
r 3 = (r°,0,0,f 3 ), n = (r?,O,fisiii0,fiCOS0), r 2 = p K ~ n - r 3 , (72) 

and R z ( y )(5) are rotations around axis z (y) with Euler angle 5. Using the Particle Data 
Group (PDG) parametrization for the phase space integrals [25], the LO contribution to 
Imv4 00+ and Im£> 00+ of the 3-pion cut diagrams so read 

7&n,LO, s 1 f , n , n , , n , r A V 2 )i \ a(p 2 ) 



lmA 00 ' + (si,s 2 ,s 3 ) = 16( - 27r )5 / dg 3 d ^2 da d cos/3 d7 gi)^,o(gi»Pi) 
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v±lmB 00 ' + (si,s 2 ,s 3 ) = 16 ( 27r )5 J dq ° d « dcos/j d7 ^4^ + J _(p x , ;±(qi,p. 



(73) 



where Ag^^g^pj), A^. are the tree level amplitudes for 7r°7r°7r + — > 7r°7r°7r + and 
7T + 7r + 7r _ — > 7r°7r°7r + . At this order the are no contribution of these topologies to Rei?oo+- 
We have computed numerically the integrals in f!73|) . The correction induced by (1731 to 
the other LO terms is always much below the 0.5%. We find so perfectly consistent to 
omit these corrections at this level of precision. Similar expressions and conclusions hold 
for the FSI in the decay of K L — ► 7r°7r 7r°. For this case one just has 

Im A)00 (Sl,8 2 ,8 3 ) = 16 ( 27r )5 J ^ dcOS ^ d ^ [^000 (Px, gi)^S,o(gi»Pi 



7 37r,LO , » 1 



u ± ImSooo"(si,s 2 ,s 3 ) = 16 (^)5 / dg " d ^2 da dcos/3 d7 A^^, g<)^3?,±(?i>Pi) 



(74) 



where A^o(?i)Pi)> ^3?r ±(?i> Pi) are ^e tree level amplitudes for 7r 7r°7r° — > 7r 7r 7r° and 



7T + 7T 7T° — > 7r 7T 7r . 
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